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Abstract. In LiHO2.Yi_2.F4, the magnetic Holmium Ho 3+ ions behave as effective 
Ising spins that can point parallel or antiparallel to the crystalline c-axis. The 
predominant inter-Holmium interaction is dipolar, while the Y 3+ ions are non- 
magnetic. The application of a magnetic field B x transverse to the c-axis Ising direction 
leads to quantum spin-flip fluctuations, making this material a rare physical realization 
of the celebrated transverse field Ising model. The problems of classical and transverse- 
ficld-induced quantum phase transitions in LiHo2;Yi_2;F4 in the dipolar ferromagnetic 
(x = 1), diluted ferromagnetic (0.25 < x < 1) and highly diluted x < 0.25 dipolar 
spin glass regimes have attracted much experimental and theoretical interest over the 
past twenty-five years. Two questions have received particular attention: (i) is there 
an antiglass (quantum disordered) phase at low Ho 3+ concentration and (ii) what is 
the mechanism responsible for the fast B x -induced destruction of the ferromagnetic 
(0.25 < x < 1) and spin glass (x < 0.25) phases? This paper reviews some of the 
recent theoretical and experimental progress in our understanding of the collective 
phenomena at play in LiHo2;Yi_2;F4, in both zero and nonzero B x . 
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Figure 1. The left panel shows the crystalline tetragonal structure of LiHoF4 
with lattice spacing a = 5.175 A and c = 10.75 A. The thin solid lines illustrate 
the supercxchangc pathways between magnetic Ho 3+ ions (black circles) mediated 
via fluorine F~ ions (green circles). In LiHo x Yi_ I F4, the magnetic Ho 3+ ions are 
substituted randomly by non-magnetic Y 3+ ions, with the lattice structure remaining 
the same for all x. The global Ising direction is along the c axis and a magnetic 
field B x is applied perpendicular to that axis in transverse field experiments. The 
right panel shows a schematic temperature (T) — dilution (x) — transverse field (B x ) 
phase diagram of LiHo 2: Yi_ 2: F4. At high temperature T or large B x , the system is 
in a paramagnetic (PM) phase. For x = 1, LHI0F4 is an Ising dipolar ferromagnet 
(FM) below T c = 1.53 K. Upon the substitution of Ho 3 + by Y 3+ , the FM phase 
persists down to x ~ 0.25. Upon cooling the diluted FM for 0.25 < x < 0.5, magnetic 
susceptibility measurements find the development of a "fcrroglass" regime (FM or 
SG?). As x decreases, random frustration builds in and, for x < 0.25, a dipolar Ising 
spin glass (SG) phase develops (shaded blue region). Whether the dipolar SG exists 
down to x = + or an exotic quantum disordered "antiglass" (AG) phase occurs at 
x > (red line segment) is not yet resolved and is the subject of much controversy. 
The application of a transverse field B x introduces quantum fluctuations causing a 
reduction of T c (x) (surface delineated by dotted lines) which, ultimately, drives a zero 
temperature quantum phase transition at zero temperature (dash-line curve in the 
x — B x plane) . 

1. The pure LiHoF 4 material 

1.1. Dipolar ferromagnetism 

The insulating rare-earth compound LiHoF 4 has the scheelite structure depicted in Fig. 
1 PQ . In this material, only the Holmium ions, Ho 3+ , are magnetic. LiHoF4 forms 
a tetragonal structure (space group C| h — IAi/a) with lattice constants a = 5.175 
A and c = 10.75 A. There are 4 magnetic Ho 3+ ions per unit cell with fractional 
coordinates (0, 0, =), (0, |, f ), (s, 5, 0) and (|, 0, |). In Fig. 1, the lines connecting the 
central Ho 3+ ion with neighbouring F~ ions and then to the nearest Ho 3+ ions illustrate 
the superexchange pathways between nearest-neighbour Ho 3+ cations [2]. As in most 
magnetic rare-earth compounds, the unfilled 4/ electronic orbitals of the Ho 3+ ions 
are not very spatially extended. This makes the exchange interaction weak and causes 
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magnetostatic dipolar couplings to be the strongest inter-ionic interactions. Indeed, 
LiHoF 4 is a dipolar Ising ferromagnet with a magnetic moment of ~ 7fiB per Ho 3+ 
ion j2] and with a critical temperature T c m 1.53 K [2]. In the mid-seventies, crystal- 
field calculations [3], susceptibility [I] and spectroscopic measurements [5] showed that 
LiHoF4 is an Ising-like dipolar ferromagnet with a highly anisotropic (/-tensor (g± « 
and g\\ ~ 14). The dipolar interaction has a long and interesting history and this physical 
realization rekindled interest in dipolar ferromagnetism. Luttinger and Tisza had found 
that the ground state of a dipolar ferromagnet depends on both the lattice structure 
and sample shape [6]. However, taking the demagnetization field into account, Griffiths 
showed that, in the absence of an applied magnetic field, the free energy is independent of 
the sample shape [7]. Using renormalization group arguments, Larkin and Khmelnitskii 
[8] and Aharony [9] showed that d* = 3 is the upper critical dimension for mean- 
field critical behaviour in a dipolar ferromagnet. This observation led to significant 
experimental efforts to find the expected logarithmic corrections to mean-field theory 
in LiHoF4 [TQl EH EE21 [131 [H] . The demagnetization field also renders the ferromagnetic 
transition quite peculiar in the presence of dipolar interactions |T| as these oppose a 
uniform magnetization and leads to the formation of domains below the transition [5lfT5]. 
The size and shape of the domains depend on an energy balance between surface and 
bulk contributions [16j. For LiHoF 4 , the ground state domain configuration is expected 
to consist of parallel antiferromagnetically aligned sheets [El [TBI US]- The size of the 
domains increases with decreasing demagnetization factor and, in the limit of infinitely 
long thin samples, a state of uniform magnetization is possible. Also, for thick samples, 
a branched pattern is predicted in the domain structure for a dipolar Ising ferromagnet 
[TT1 [TH] . LiHoF4 is an excellent test case for theories of domain formation due to the 
availability of an accurate microscopic model of the material (see below). It remains an 
open experimental challenge to verify the nature of the predicted rich domain structure. 

1.2. Microscopic model and transverse field Ising model for LiHoF^ 

The magnetic properties of LiHoF4 are determined by the 4/ 10 electrons of the Ho 3+ 
ions. Applying Hund's rules results in a 17- fold degenerate 5 /g electronic ground state. 
The crystal field from the surrounding ions, described by a crystal field potential Vcf, 
lifts this degeneracy and gives rise to a two-fold non-Kramers degenerate ground state, 
and an excited singlet, \ipe), at an energy of approximately 10 K above \?Po)- 
The minimal microscopic Hamiltonian in the presence of a transverse magnetic field 
B = B x x can be written as |20j : 



where fj,,u = x,y,z. C%j is the magnetic dipole interaction, = [S^ u \fij\ 2 — 

3(r*jj) M (f*jj) iy ]/|fi : ,| 5 , J ex is the nearest-neighbour exchange interaction. A is the strength 
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of the hyperfine interaction and Ii (I = 7/2) is the total angular angular momentum 
vector of the Ho-nucleus at the i-th site. Disregarding the hyperfine interactions, an 
effective spin 1/2 model for LiHoF4 from the above microscopic model of Eq. (TjQ) was 
derived in Ref. [20] and revisited in Ref. [21]. Without hyperfine interactions, the 
single site Hamiltonian (first two terms in Eq. (CQ)) can be diagonalized numerically for 
arbitrary value of the transverse field B x . The ground state doublet is split by B x , 
with an energy A(B X ) between the two states. Meanwhile, the excited \ip e ) singlet state 
remains more than 10 K above the split doublet. Consequently, the Jj operators can 
be projected onto the two-dimensional subspace of the two lowest energy eigenstates 
via the relationship J? = C^o + Yl v = x y,z £W(-^) a T) hence providing a description 
in terms of effective pseudospin-1/2 operators. The strongest effective interaction is 
oc (C zz ) 2 Cffa^a? and, to an accuracy of a few percent [20j [21], the effective model can 
be written (up to a field dependent constant) as 

Rising = - £ of + \{g^ B C zz f £ + J^{C zz f J2 ofoft2) 

Note that the C^ v parameters and the effective transverse field r, T = A/2, depend on 
the applied physical transverse field B x . In the limit of small B x , A(B X ) oc B 2 X [201 12T] . 
In Eq. ([2]), A/2 plays the role of an effective transverse field acting on the x-component 
of the pseudospin <?, — hence the realization of a transverse field Ising model (TFIM) 
in LiHoF 4 [22]. 

For B x = (A = 0), the effective dipolar model (T2]) has been analyzed using 
mean-field theory j2j H] and classical Monte Carlo simulations. Numerical Monte 
Carlo simulations are complicated by the difficult nature of the conditionally-convergent 
dipolar lattice sums. An early Monte Carlo simulation using free boundary conditions 
for rather small system sizes found a critical temperature of 1.89 K [23]. To obtain 
better accuracy with long range dipolar interactions, there are two standard approaches 
that can be implemented: Ewald summation [191 Ell [Ml [251 EH [27] and cavity methods 
[201 [28] . An early Monte Carlo simulation using the Ewald summation method to 
a dipolar Ising ferromagnet on a body-centered lattice was carried out by Xu et al. 
[2H[25]. Monte Carlo simulations using the cavity method [20] and the Ewald method 
[T9l I2T] have determined the critical temperature for the dipolar model in Eq. ([2]) with 
A = on the LiHoF 4 lattice and estimated the nearest-neighbour (antiferromagnetic) 
exchange interaction J ex needed to lower the T c Monte Carlo value in order to match it to 
the experimental value of 1.53 K [2]. A more recent study, using the Ewald summation 
method and system sizes up to 32 000 dipoles, found that the magnetization, specific 
heat and susceptibility do match experimental results at a quantitative level [19]. The 
pure LiHoF 4 material can therefore be viewed as one of the best realizations of a dipolar 
Ising model. A review of materials that behave like Ising systems can be found in 
Ref. [29]. 
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1.3. Hyperfine interaction 

As in other Ho-based materials [30j [3T] , there is a strong hyperfine coupling (A = 0.039 
K) between the electronic (J = 8) and nuclear (7 = 7/2) moments in LiHoF4. The 
hyperfine contribution dominates the specific heat below 0.5 K [2], and leads to a 
substantial increase in the critical transverse field at low temperatures [22]. As we shall 
discuss in Section [21 the hyperfine interaction also strongly influences the phase diagram 
of the dilute compound LiHo x Yi_ x F4 in a transverse field [321 [33] . The hyperfine 
interaction is therefore important and may not be omitted in a proper quantitative 
description of the low-temperature properties of LiHo x Y 1 _ x F 4 . 

1-4- Transverse field effects 

At a continuous phase transition, the order parameter fluctuates coherently over 
increasing length scales as the critical point is approached. At finite temperatures, 
the fluctuations are classical in nature sufficiently close to the critical temperature. In 
some systems, at temperatures near absolute zero, a phase transition can be induced 
by tuning an external parameter such as a magnetic field, pressure, doping or disorder. 
In such a case, the fluctuations are quantum mechanical in nature and the transition 
is called a quantum phase transition (QPT) [311 EH]. The archetypical model for a 
quantum phase transition is the transverse field Ising model (TFIM) [36|l37]. There are, 
however, not many real Ising magnets with experimentally accessible QPTs [38] since 
the required transverse field usually far exceeds laboratory fields (currently less than 50 
T for DC fields). This would be the case for dipolar Ising spin ice materials [39], but 
not for the Dy(OH)3 and Ho(OH)3 dipolar Ising ferromagnets [38]. Due (i) to the weak 
dipolar interaction and (ii) the low- lying singlet |^ e ) above the ground Ising doublet 
|^), LiHoF4 is one of very few magnetic systems where one can observe the vanishing 
of the magnetization as a transverse field is applied [22] . Interestingly, recent neutron 
scattering measurements of LiHoF 4 have revealed that the relatively strong hyperfine 
interactions in Ho mask the true QPT of the hybridized electronuclear critical mode, 
and the energy gap in the solely electronic degrees of freedom remains finite as the 
system passes through the QPT [30] . 

Given the accuracy of the effective Ising model for LiHoF 4 in the absence of a 
transverse field, the material would seem to constitute an ideal testing ground for 
the inclusion of non-commuting quantum-mechanical terms in realistic spin models of 
magnetic materials. As outlined above, there is a well defined procedure to describe the 
effect of the transverse field B x in LiHoF4 within a TFIM [201 El] - However, a quantum 
Monte Carlo calculation [20] of the critical temperature, T c , vs B x phase diagram reveals 
only a qualitative agreement with the experiment of Ref . [22] . In the low temperature 
region (below 0.5 K), the effect of the hyperfine interactions on the B x —T phase diagram 
becomes important, and this is only approximately taken into account in the simulation 
of Ref. [20] . Recent theoretical work proposes an alternative way to incorporate the 
effect of the hyperfine interactions in an effective model of LiHoF 4 in nonzero B x [33] . 
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Notwithstanding the importance of hyperfine effects at low temperature, one notes that 
even close to the classical critical temperature of T c = 1.53 K, the experimental [22] 
phase boundary rises much more steeply than in the model of Eq. (T2]) [201 EI]- To 
compound the puzzle, a very recent experiment using dilatometry [4T] confirms the 
previous T C (B X ) experimental phase diagram obtained on the basis of susceptibility 
measurements [22]. The reason for this stability of the ordered state with respect to 
weak B x is an open question. Possible explanations include magnetostriction effects, 
small terms that are neglected in the effective model of Eq. (j2J) (see Refs. j2(Jl ETJ 
for a discussion of those terms) or anisotropic exchange terms, such as ^(r^J^Jj, 
altogether ignored in the microscopic model of Eq. (CO). While the latter interactions 
may largely be inconsequential at the classical transition, when projected within the 
manifold of |^) doublets when B x = [42], these could perhaps end up affecting the 
T C (B X 7^ 0) phase diagram when B x ^ 0. 

2. LiHo a ,Y 1 _ a .F 4 — transverse field and random field physics 

2.1. Disordered magnetic systems 

All real materials contain a certain amount of frozen-in random disorder such 
as impurities, interstitials and crystalline defects. One may then ask how the 
thermodynamic properties and the stability of the ground state of an otherwise idealized 
pure system are affected by various types of random disorder [431 |44]. Since the mid 
1970s, magnetic systems have proven to be almost ideal paradigms to study the role 
of random disorder in condensed matter systems. Examples of cornerstone questions 
in the field of random disordered systems are: Does a transition that is first order 
in a pure system remains first order in the presence of random disorder [IS]? For a 
second order transition, does an arbitrary small amount of random disorder change the 
critical exponents and hence modify the universality class [46]? Does the long range 
order exhibited by a pure system survive arbitrary small random fields conjugate to the 
order parameter [47]? When the disorder is large enough and the interactions are in 
strong competition and lead to random frustration, is there a thermodynamic transition 
at nonzero temperature to a spin glass state [HI 09]? Under what conditions does a 
random magnet exhibit nonanalytical behavior above its critical temperature, in the 
so-called Griffiths phase, which is a remnant of the long range order phase displayed by 
the parent disorder-free system [50]? The above questions were originally by and large 
investigated for classical phase transitions at nonzero temperature. In the late 1980s, 
those and similar questions started to attract considerable interest in the context of 
quantum phase transitions [35]. Here again, the LiHoF4 material, or more precisely, its 
Y 3+ diamagnetically diluted LiHo x Y 1 _ 2 .F 4 variant, has proven to be a highly interesting 
material to investigate quantum phase transitions in a disordered systemt and tuned 
by an external parameter. As in the pure LiHoF4 material, the parameter tuning the 
level of quantum (spin flip) fluctuations in LiHo x Y 1 _ ;l E4 is a magnetic field B x applied 
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perpendicularly to the Ho 3+ Ising spin direction. 

2.2. Transverse field and random field effects in LiHo x Yi_ x F^ 

The LiHo x Yi_2;F4 compound forms a solid solution from the pure dipolar ferromagnetic 
LiHoF 4 to the non-magnetic LiYF 4 without a change of the crystalline structure and, 
thanks to the close ionic radius of Ho 3+ and Y 3+ , LiHo x Yi_ :r F4 displays a minimal 
dependence of the a and c lattice parameters on x. The critical ferromagnetic transition 
temperature, T c (x), decreases as the magnetic Ho 3+ ions are substituted by non- 
magnetic Y 3+ . For 0.25 < x < 1, LiHo z Yi_ x F4 can be described as a diluted ferromagnet 
(see Fig. 1). As discussed in Section [LTl the main interactions between the Ho 3+ ions are 
magnetostatic dipole-dipole interactions. The dipolar interaction is inherently frustrated 
since its sign, that of in Eq. ([I]), depends on the orientation of the magnetic moments 
with respect to = fj — the position vector of ion j relative to ion i. Hence, as 
a result of the random substitution Ho 3+ — > Y 3+ , random frustration builds up as x 
decreases until, for x ~ 0.25, the diluted ferromagnetic state gives way to a dipolar spin 
glass phase (see. Fig. 1). We note that there has been much debate as to whether such a 
spin glass phase occurs for x < 0.25 in LiHo rr Y 1 _ a ,F 4 [51], E2]. We proceed in this section 
with the assumption that a spin glass phase does exist, at least for 0.15 < x < 0.25. We 
return in Section [3] to the question of the existence of a thermodynamic spin glass phase 
in LiHo x Y 1 _ ;c F4 for x < 0.25, below the concentration where the diluted ferromagnetic 
phase has disappeared (see Fig. 1). 

The essence of the research on transverse field induced quantum fluctuations and 
quantum phase transitions in LiHo a; Yi_ :r F4 originally lay in the expectation that this 
material may be a physical realization of a transverse field Ising model (TFIM) with 
Hamiltonian 



with effective transverse field T. In Eq. ([3]), the randomly frustrated fl^ couplings 
model, or rather mimic, the site-diluted £f? — > €i€jC*j in Eq. (T5]), where e% = \ when 
the site i is occupied by a magnetic Ho 3+ ion and = when it is occupied by non- 
magnetic Y 3+ . In theoretical studies, important simplifications have been made by 
taking a Gaussian distribution for Q^J of mean f2o and standard deviation Q and setting 
€i = 1 Wi [531 EH ESI EH]- Of particular interest in these studies was the prediction 
of Griffiths phase physics near the QPT both in two and three dimensional spin glass 
versions (f2 = 0) of Eq. ([3]) [55] . 

Two puzzles were identified when the behaviour of LiHo x Yi_ x F4 subject to B x ^ 
was first compared with expectations from theoretical and numerical studies. Firstly, 
while dT c /dx oc x in the diluted ferromagnetic (FM) regime (0.25 < x < 1.0, see 
Fig. 1), as in mean-field theory [57] . the rate at which T C (B X ) is reduced by B x 
gets progressively faster than mean-field theory predicts as the Ho 3+ concentration x 
decreases [58]. This means that, in comparison with the energy scale for FM order 





(3) 
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that determines the classical transition at T C (B X = 0,x), B x becomes more efficient 
at destroying FM order the lower the concentration x is [58]. Secondly, for B x = 0, 
LiHoo.i67Yo.833F4 has generally been believed to exhibit a conventional spin glass (SG) 
transition [51] (see however Ref. [52j [59]), with the transition signaled by a nonlinear 
magnetic susceptibility, x 3 , diverging at T g as XsiT) oc (T — T 9 )~ 7 [48J. However, 
X 3 (T) in LiHo .i67Y .833F4 becomes less singular as B x is increased from B x = 0. This 
suggests that, in contrast with numerical work [54| [55] , no ^-driven quantum phase 
transition between a polarized quantum paramagnetic state and a SG state occurs as 
T [601 EH- 

The resolution of these two puzzles lies in the observation that, in the properly 
derived effective low-energy theory akin to Eq. ([2]) when the occupation factor q ^ 1 
[62} [63], the transverse field B x induces not only quantum fluctuations, but also random 
fields [32], [HJ, hi, that couple linearly to the z component of <?.; [321 [331 E21 E31 EH 165] . 
There are two key physical ingredients responsible for these random fields. Firstly, the 
magnetic field B x , perpendicular to the c-axis Ising direction, induces a finite expectation 
value of the Ho 3+ magnetic moments along x, (/i^) = g/iB(Jf), at all temperatures. 
Quantum mechanically, this moment arises from the S x -induced admixing of the 
(initially strictly Ising) |^) ground doublet with the excited crystal field states. 
Secondly, the bare microscopic Hamiltonian in Eq. ([1]), as opposed to the effective 
TFIM of Eq. ([2]) and its simplified version in Eq. ([3]), contains off-diagonal (anisotropic) 
C!£j (fi ^ v) terms of dipolar origin that couple the various components Jf of Jj. 
In pure LiHoF4, perfect mirror symmetries make the lattice sum h\ = Y2j (Jj) 
vanish identically. For diluted LiHo x Y 1 _ a .F 4 , h\ no longer vanishes; the transverse 
magnetic moments (/xf) induce, via the random ej>CfJ{Jf ) terms, internal random fields 
that act on the z component of Jj, Jf . In other words, diluted LiHo x Yi_ x F 4 subject 
to a transverse field B x maps onto a TFIM with the addition of correlated random 
fields [621 E3]. For 0.25 < x < 1, this makes LiHo 3; Y 1 _ :c F4 a unique example of a 
random field Ising ferromagnet model. Indeed, the experimental procedure for creating 
such a system is convoluted and typically considers a diluted Ising antiferromagnet in 
a uniform longitudinal field B z [661 EZ] • Analytical [331 E21 E31 EU [65] and numerical 
[621 163"! E3J studies have shown, rather convincingly, that the presence of -B x -induced 
random fields h\ can explain the two aforementioned puzzles observed in LiHo x Yi_ a; F4 
in nonzero B x , in both the diluted FM regime and in the SG regime. In particular, the 
-B^-induced random fields hi eliminate the divergence of the nonlinear susceptibility % 3 
in the spin glass regime x < 0.25 [321 E21 164] . 

2.3. Open questions regarding the random field physics in LiHo x Yi_ x F± 

Because of space constraint, we have focused above on the role of the induced random 
fields on the phase transitions that occur in the diluted ferromagnetic and spin glass 
regimes. As the critical (ferromagnetic or spin glass) temperature decreases, and drops 
below approximately 0.5 K, the role of the hyperfine interaction in inhibiting quantum 
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fluctuations and strongly renormalizing the critical transverse field cannot be ignored 
[32| [33] . While progress has been achieved as to how to proceed [33] , work does remain 
to be done to investigate via quantitatively accurate calculations (e.g. quantum Monte 
Carlo simulations) the role of hyperfine interactions on the global x — T — B x phase 
diagram of LiHo x Yi_a;F4. 

In the context of the physics of induced random fields, it is worthwhile to note 
that theoretical work predicts that in the simplest TFIM of Eq. ([3]), the addition of 
uncorrelated random fields should destroy the zero temperature quantum criticality 
and drive the renormalization group flow towards the fluctuationless (classical) zero 
temperature random field fixed point [68]. It would be interesting to investigate the 
implication of that result for the correlated h\ random fields discussed above and that 
exist in LiHo 2 ,Yi_ x F4 in nonzero B x . Finally, we note that the presence of induced 
random fields has been invoked to interpret the transverse field and temperature 
dependence (e.g. scaling behavior) of the magnetic susceptibility in LiHo x Y 1 _ :c F4 
(x = 0.44). In particular, these random fields have been postulated to give rise to 
Griffiths' singularity effects in the paramagnetic regime at temperatures just above 
T c [69]. It would be useful to explore this claim further. In fact, a systematic 
experimental investigation of the fascinating physics of random field systems in the 
diluted ferromagnetic regime of LiHo x Y 1 _ a .F4 is probably warranted. 

3. LiHo x Yi_ x F 4 — spin glass vs antiglass 

3.1. Spin glass vs antiglass physics in LiHo x Y\_ x F^ 

The previous section discussed the properties of LiHo x Y 1 _ x F4 in nonzero transverse field 
B x . We now turn to the question of the existence of a spin glass phase in LiHo x Yi_ 3; F4 
for x < 0.25 and in zero B x . 

The x — T phase diagram of LiHoaYi_ x F4 in the limit of large dilution (x < 0.25) 
is a subject of much controversy. Theoretical studies of the dilute Ising dipolar model 
predict that a spin-glass state may be favored over ferromagnetic order at high dilution 
[241 170] . Since there is no percolation threshold for long range dipolar couplings, a 
thermodynamic spin glass state is expected to persist all the way to x = + [70J. 
Furthermore, a study of the quantum corrections to the Ising model arising from 
interaction-induced virtual crystal field fluctuations finds that these quantum effects 
should be small in LiHo x Y 1 _ a .F 4 [42]. For example, this is in contrast to the Tb 2 Ti 2 07 
and Tb2Sn20"7 pyrochlore Ising magnets [TTJ [72] and, therefore, are unlikely to destroy 
the spin glass phase in LiHo x Y 1 _ a F 4 at finite Ho 3+ concentration. Consequently, in 
light of the accumulated evidence for a thermodynamic spin glass phase transition in 
the three-dimensional Edwards- Anderson (EA) Ising spin glass model [49J, one would 
also expect such a transition in LiHo x Yi_ x F4 for < x < 1 since the EA model and 
dipolar systems should be in the same universality class for a spin glass transition [73] . 

Measurements of the imaginary part of the magnetic susceptibility, x\ m 
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LiHo .o45Yo.955F 4 by Reich et al. [73], [75] was the first to challenge the above theoretical 
picture of a spin glass phase down to x = + . In a typical glassy system, the 
frequency distribution of \ broadens upon lowering the temperature [IB]. In contrast, 
the authors of Ref. [TjJ [75] found a distribution that narrows as the temperature 
is reduced, indicative of an overall speed-up of the dynamics and an "antiglass" 
(quantum disordered) ground state. The measurement was later repeated [76] to 
lower temperatures (50 mK), still showing an evident narrowing of the spectrum. 
The manifestation of dissipationless dynamics of a broad inhomogeneous distribution 
of quantum oscillators within the inferred antiglass regime was further demonstrated 
through (i) hole burning in the AC susceptibility spectrum [76] and, (ii) perhaps most 
interestingly, via the observation of persistent oscillations for up to 30 seconds after an 
external AC field is turned off [76]. It has been argued that this antiglass behaviour may 
originate from a quantum mechanical entanglement of the magnetic dipole moments and 
the formation of a sort of random-singlet state [77]. However, as mentioned above, the 
quantum effect corrections to the otherwise perfect Ising nature of the single-ion crystal 
field ground doublet is expected to be very small in LiHo x Y 1 _ a ,F 4 [32]. Also, while 
hyperfine interactions do strongly renormalize the critical transverse field in transverse 
field experiments (see Section 2) [32j [33], at the same time, preliminary calculations 
suggest that these interactions may not be efficient at quenching the spin glass phase at 
a Ho 3+ concentration of order of 5% [32] • It is therefore not clear how the entanglement 
mechanism proposed in Ref. [77] on the basis of an ad- hoc model can really proceed 
efficiently when starting from the microscopic model of Eq. (JTJ). 

Apart from the question of the existence of an antiglass phase, there have also 
been contradictory experimental reports on the spin glass transition in the regime 
0-1 % x ~ 0-2, where it was generally believed to be well established. Specifically, 
Jonsson et al. found that the non-linear susceptibility % 3 does not diverge for samples 
with x = 4.5% and, surprisingly, for x = 16.5% as well. These results prompted 
Jonsson et al. to conclude that there is no thermodynamic spin glass phase at all in the 
compound [52]. In response to this work, Ancona- Torres et al. re-examined the non- 
linear susceptibility for x = 16.7% and x = 19.8% and concluded that X3 does diverge 
|51j . A possible explanation put forward in Ref. [51] for the discrepancy between the 
works of Ref. [52] and Ref. [51] is that the sweep rates of the magnetic field in the work 
of Jonsson et al. [52] were too fast and the fields too strong to expose the true spin 
glass critical behaviour of the system. To compound the debate, Jonsson and co-workers 
[59] continue to maintain that evidence for a spin glass transition in LiHo a; Y 1 _ :r F4 for 
x = 16.5% is inconclusive, unlike the claims of Ancona- Torres et al. [5Tj. 

Several numerical studies have aimed at finding a spin-glass transition in effective 
dipolar Ising models mimicking the dilute compound LiHo x Y!_ x F4. The original 
work of Xu et al. only identified the disappearance of the long range ferromagnetic 
order for a diluted dipolar body-centered cubic lattice and did not investigate the 
spin glass freezing [23]. The traditional way to identify the glass transition has 
been to calculate the Edwards-Anderson overlap between two replicas, #1 and #2, 
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Figure 2. The left panel shows a plot of the imaginary part of the AC susceptibility, 
X" , as a function of frequency, /, normalized to its maximum value ^" max 

at frequency 

/max, and for temperatures between 77 mK and 260 mK. The figure shows a broadening 
of the spectrum below and above the peak frequency / ma x at which x"(/) peaks 
(adapted from Ref. [79] )• The right panel shows the real part of the AC susceptibility, 
%'(/), as a function of temperature. One notes the general good agreement with the 
lim/^ox"(/) data set from Ref. [75] down to a temperature ~ 180 mk as well as 
with the DC susceptibility data determined via classical Monte Carlo [27]. The figure, 
adapted from Ref. [79] . illustrates the disagreement with the results of Jonsson et al. 
in Ref. [52] as well as with the T -0 75 behaviour reported in Ref. [77] and argued to 
support an entangled collective singlet ground state. 

q = 1/N of 1^25 an d to determine the temperature at which the corresponding 
Binder ratios, g — 1 — [(<7 4 )]/3[(g 2 )] 2 , for different system sizes intersect. However, 
several Monte Carlo studies [26] 127] have been unable to detect a spin glass transition on 
the basis of such intersection, or other finite-size scaling properties of q [78]. These "first 
generation" Monte Carlo simulations of diluted dipole-coupled Ising spins did therefore 
suggest that the absence of a thermodynamic spin glass transition in LiHo x Yi_a;F4 at 
small x could possibly have a classical origin. However, as we discuss in the following 
section, more recent experimental and numerical studies reveal a different picture. 

3.2. Recent developments in spin glass physics in LiHo x Yi_ x F4 

Since much of the evidence for the antiglass phase is based on the scaling of the 
imaginary part of the AC susceptibility, the reproducibility of these measurements is 
very important. Recent measurements by Quilliam et al. [79] contradict the earlier 
results [7_U ESI EE] and find that the frequency distribution of x'\ broadens as the 
temperature is reduced also in the high dilution limit (x = 4.5%) as shown in the left 
panel of Fig. 2. Furthermore, the quasi-static susceptibility, Xo ; reported in Ref. [79] 
does not show the unusual T -0 75 divergence found earlier [77]. We note, perplexingly, 
as noted in Refs. [27J ES], that xo reported by Ghosh et al. in Ref. [77] does not agree 
with the previous results [74] EH] by the same group on, presumably, the very same 
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sample. Meanwhile, the experimental Xo data of Quilliam et al. [79] agree reasonably 
well with the early results of Refs. [7U [75] down to a temperature of approximately 180 
niK, as well as with the DC susceptibility found in a classical Monte Carlo simulation 
of Eq. (j2J) with B x = [27], as shown in the right panel of Fig. 2. The work of Ref. [79] 
would appear to give a strong indication that there is conventional spin glass freezing 
in the high dilution limit of LiH0a.Yi_2.F4. Furthermore, the specific heat data reported 
in Ref. [80] do not display any of the unusual features that were previously interpreted 
as a telltale signature of entangled magnetic moments within a quantum disordered 
antiglass state [77] . Note that the specific heat determined via classical Monte Carlo 
simulations [27] does not agree particularly well with either set of experimental data. 
This discrepancy may be in part due to the difficulty in taking into account accurately 
of the hyperfine contribution to the experimental specific heat data. 

The reason for the difference between the AC susceptibility and specific heat 
experimental results of Refs. [7JJ [75j CEJ [77] and Refs. [791 EQ] remains an open question. 
Due to the very slow dynamics, all these measurements are technically very demanding. 
However, we believe that it would be important to settle the question of the nature of 
the high dilution phase in LiHozYi-aT^ Perhaps the study of other highly diluted Ising 
magnetic materials could help shed some light on the physics at play in LiHo x Y 1 _ :r F4. In 
this context, a study of the highly diluted Ising spin ice materials (Dy x Y 1 _ x ) 2 Ti 2 07 and 
(Ho x Yi_ x ) 2 Ti 2 7 [M] as well as the proposed Ho x Yi_ x (OH) 3 and Dy x Yi_ x (OH) 3 TFIM 
materials [38] to compare with LiHo x Y!_ x F 4 could prove interesting and instructive. 

The failure of previous Monte Carlo simulations [26j[27J,[78] to find a thermodynamic 
spin glass transition in the dipolar model would seem to have been explained recently. 
As in three-dimensional Edward- Anderson Ising spin glass model [49], the spin glass 
correlation length £ sg turns out to be a much better indicator of a spin glass transition 
than the Binder ratio g. A recent Monte Carlo study [82] that focuses on £ sg confirms 
that this is also the case also for the dipolar Ising model. In Ref. [82], a finite-size 
crossing of £ sg divided by the linear system size L, £ gg /L, is observed, providing good 
evidence for a finite temperature spin-glass transition, in agreement with the most recent 
AC susceptibility experiments [79] . However, the critical temperature determined by the 
crossing of £, S g/ L is slightly different for £ sg along the a and c axes. This indicates that 
the finite size corrections to scaling are not under control for the system sizes considered 
and simulations on larger system sizes are therefore needed. A very recent numerical 
study argues for a quasi-long-range ordered spin glass phase in a system of diluted Ising 
dipoles [83J. 

3.3. Open questions regarding the spin glass physics in LiE.o x Y\_ x F^ 

The single most important open question is the nature of the magnetic state in the 
limit of high dilution (x < 0.1). Is it a traditional spin glass, as suggested by the 
recent measurements of Quilliam et al. [79] or is it an unusual antiglass spin liquid 
[74^ [75] [76] [77] ? One may even ask whether there is a spin glass phase at all for 
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x < 0.25 [5TJ E2J ES]- It would appear that the only way to resolve this issue is through 
more experiments. The "sample quality" (e.g. impurity, homogeneous Y 3+ disorder, 
etc) may be of importance in this highly dilute limit and this issue should be examined 
carefully. Sample shape and demagnetization effects should not affect the existence of 
the spin glass phase, but this should nevertheless be investigated. There is a significant 
qualitative difference between the static and dynamic susceptibility measurements 
reported by three research groups, and a new and independent series of measurements 
could help shed light on this outstanding question. It appears that the nonlinear 
susceptibility results of Ref. [52] are being dismissed on the basis of measurements 
that do not really probe the spin glass critical regime. However, disagreement on this 
interpretation remains [59]. The hole-burning and persistent oscillations observed by 
Ghosh et al. [76] are also indicative of a very unusual physical state, yet they have not 
been reported or reproduced by another group. To reiterate, we believe that at this 
stage, in order for the field to move forward, more experiments that carefully explore 
the critical regime, notwithstanding the highly technical burden associated with those, 
are very much needed. 

On the numerical side, it would be of interest to go beyond the classical Monte 
Carlo simulations that have so far been applied to the dilute system. The effects of 
the hyperfine interactions and off-diagonal interactions between the pseudospins are 
challenging to incorporate in a quantum Monte Carlo simulation since they lead to a 
sign problem. Nevertheless, it would be important to invest some effort in this problem 
in order to achieve more quantitative comparisons with experiments. 

4. Conclusion 

Due to the existence of a well-established two-state Ising model for the rare-earth 
compound LiHo x Y 1 _ a .F4, this material constitutes a model magnet where theory, 
experiments and simulations can be compared and contrasted at a quantitative 
level. Despite the apparent simplicity of the underlying model for the parent 
compound, the introduction of disorder through dilution and quantum effects through 
a transverse magnetic field leads to a complex system which continues to challenge 
experimentalists and theorists alike. The two foremost outstanding questions are 
whether a thermodynamic spin glass phase survives in the limit of high dilution of 
the magnetic Ho 3+ ions and the reason for the unexpected stability of the ferromagnetic 
ground state of pure LiHoF4 with respect to an applied transverse field. Continued 
research into this remarkable realization of the Ising model will hopefully continue to 
shed light on such diverse topics as random fields, coherent oscillations and tunable 
quantum fluctuations. 
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